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ABSTRACT 


Bourdon tubes are the most commonly used elastic elements in mechanical 
pressure gauges. Basic theory describing the principles of Bourdon behavior is readily 
available but very few analytical studies have been published that model Bourdon 
element behavior. The purpose of this research is to develop an analytical model to 
determine tip displacement in tubular elastic elements using basic principles of solid 
mechanics. To determine the validity of the results obtained from the analytical model, a 


study of results from finite element solutions and experimental data is also presented. 
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I, INTRODUCTION 


A. MOTIVATIONS FOR THIS RESEARCH 

Bourdon tubes, the elastic elements used in mechanical pressure gauges, are 
quite popular, even today, due to their simplicity. When a Bourdon tube is fabricated 
from conventional alloys and operated within its elastic limit, the tip of the tube moves a 
small distance. However, if the Bourdon element were fabricated from a superelastic 
alloy capable of withstanding strains on the order of 8%, strains much higher than what 
ordinary alloys can withstand, the Bourdon tube could be used in unconventional 
applications as a robotic finger or as an actuator. With these unique applications in 
mind, the study of the elastic deformation of a Bourdon tube was initiated. 

The efforts of a literature survey, the results of which are discussed later in this 
chapter, yielded little information in terms of analytical solutions to Bourdon tube 
behavior. In industries that manufacture Boaniis tubes, tip travel is obtained from 
experimental data in the absence of a mathematical model. This further prompted the 
development of an analytical solution for tip travel in a Bourdon tube. 

With today’s micro-chip technology, many people think that the trend for pressure 
gauges is moving in the direction of the electronic transducer. Surprisingly, precision 
pressure gauges demonstrate consistent and accurate behavior using Bourdon tube 
technology. As a result, the demand for Bourdon gauges is quite significant. The most 
obvious advantages characteristic of Bourdon gauges are their (1) broad pressure range, 
(2) low cost compared to electronic gauges, (3) ease of manufacturing, (4) ease of 
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operation, (5) ease of troubleshooting, (6) good accuracy (up to 0.1 % full scale 
deflection), and (7) lack of an external power source which precludes susceptibility to 
line noise and power outages. [Ref. 1,2] 

Bourdon tubes are extremely popular as elastic elements in pressure gauges. 
Despite their popularity little analytical work has been performed or published modeling 


their behavior. The work presented in this thesis provides such an analysis. 


B. BOURDON TUBE BEHAVIOR 

1. Developmental History 

In 1849, a French engineer named Eugene Bourdon invented the C-shaped 
element commonly known today as the Bourdon tube or C spring. While witnessing the 
construction of a steam engine, Eugene Bourdon noticed that the helically wound coil 
used to condense the steam was flattened during formation. To correct the deformation, 
the tube was plugged at one end and pressurized at the other. As a result, the coil 
partially unwound as the flattened tube regained its circular cross-section. Intrigued by 
what he saw, Eugene Bourdon conducted his own experiments and eventually invented a 
pressure gauge using the tip deflection of a curved tube with an elliptical cross section as 
a pressure measuring device. [Ref. 1] Figure 1 offers an exploded view of a typical C- 
type pressure gauge. 

2. Deflection Theory 

The Bourdon tube, henceforth referred to as the Bourdon, is the actuating element 


in many pressure measuring devices. Typically, the Bourdon elastic element has an oval 
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Figure 1. C-type Bourdon gauge 




















or elliptical cross-section and is coiled in the shape of a C, a spiral, or a helix as shown in 
Figure 2. One end of the tube is sealed and the other is connected to a capillary tube 
through which pressure is applied. Any pressure within the tube in excess of the external 
or atmospheric pressure causes the flattened element to elastically change its shape to a 
more circular cross-section. The resulting change in the cross-section causes the coiled 
tube to straighten. With one end of the tube fixed, the straightening effect causes the free 
end to move a small distance. This small, unenhanced movement is called tip travel as 
shown in Figure 3. [Ref. 1,3,4,5,6] 

In a C-type Bourdon, the elastic deformation of the cross-section creates 
compressive stresses along the inner wall and tensile stresses along the outer wall as 
defined in Figure 3. When added to the tensile stresses existing in both walls due to the 
applied internal pressure, a stress imbalance is created such that the tensile stress in the 
outer wall is greater than the tensile stress in the inner wall. In the process of 
straightening, the tube experiences tensile stresses in the inner wall and compressive 
stresses in the outer wall, just the opposite caused by the change in cross-section. The 
effect of straightening relieves the unbalanced stress condition created by the change in 
cross-section. The free end of the Bourdon moves with the application of pressure until 
opposing stresses within the tube reach equilibrium. [Ref. 1,4,5] 

3. Design Criteria 

As stated earlier, the first Bourdon pressure gauge was patented more than 100 


years ago. But despite its age and the fact that “technology” is now much more 








GAUGES — PRESSURE INDICATING 
ASME 840. 1-199% ; OlAL TYPE — ELASTIC ELEMENT 





C-Type Bourdon Tube Helical Bourdon Tube 


Spieal Bourdon Tube 


Figure 2. Common Bourdon elastic elements 














| 













DIRECTION 
OF 
MOTION 
RP eieiaee 
Boga) ete | 
Fo SVS 
7 Yor 
OUTER WALL P TIP 
(TENSILE STRESS) /* See 







INNER WALL 
(COMPRESSIVE STRESS) 


(a) 


eee MAJOR AXIS 








VLLLLLLLLL 20 D>» 


Lo 

2, Sa 
j MW 2x 
Na, Vp St 


SLL LLL 
SECTION A-A | 


(b) 








Figure 3. Description of tip travel in a Bourdon element (from Ref. 1) 

















advanced, very little about the design of a Bourdon gauge has changed. The factors 
considered most important in the design of tubular elastic elements are (1) the element 
material, (2) the coil length, (3) the wall thickness of the tube, (4) the cross-sectional 
shape, and (5) the desired tip travel. 

Bourdons may be made of ferrous or nonferrous metals depending on the 
application of the gauge. Admiralty brass, phosphor bronze, and beryllium copper, which 
are used to make grade B accuracy gauges, higher quality grade A and AA gauges, and 
high pressure gauges, respectively, are examples of nonferrous metals applications. Ni- 
Span C is a popular but expensive ferrous material used for its exceptional temperature 
behavior. [Ref. 1] 

As stated earlier, Bourdon elements may be coiled in the shape of a C, a spiral, or 
a helix. The shape and coil length of the elastic element depends on the pressure range 
of operation. For pressures up to 60 bar (6 MPa, 870.24 psi), C-type elements are 
typically used and for pressures above 60 bar, spiral or helical elements are used. 

[Ref. 5] 

The geometry of the cross-section, one factor affecting tip displacement, is 
chosen depending on the pressure range of operation. The distinction is made between 
low and high pressure measuring elements for the proper choice of the wall thickness. 
Low pressure gauges are loosely defined as those operating from 12 - 500 psi and 
typically have flatter cross-sections and thinner walls than high pressure gauges which 


have rounder cross-sections and thicker walls and may read pressures as high as 100,000 














psi. [Ref. 6] Low pressure Bourdons experience greater tip travel for a given pressure 
whereas high pressure Bourdons exert higher tip forces. [Ref. 1] 

Tip travel in a Bourdon element varies widely with the range of the gauge. Low 
pressure gauges, with thinner walls and smaller minor axes, are less effective in 
overcoming the frictional forces within the amplifying movement. Higher tip travel and 
smaller amplification means smoother pointer motion but also produces excessive 
stresses in high pressure elements. As a result, high pressure gauges have thicker walls, 
larger minor axes, and high amplification movements to keep tip travel small. Gauge 
manufacturers design their Bourdon elements for specific magnitudes of tip deflection. 
For ease of manufacturing and maintenance, variations in design are kept to a minimum. 
As a result, the major distinction between gauges is found within the movement 
mechanisms. 

When designing their gauges, manufacturers either determine the appropriate 
length of the sector tail and the gear ratio necessary to produce the pointer rotation for 
the desired tip travel or determine the necessary tip travel to operate a specific 
movement. [Ref. 1, 2, 3, 4,5] Bourdon tip displacement is generally linear with pressure 
such that 25%, 50%, and 75% of the total applied pressure is equal to 25%, 50%, and 
75% of the total tip displacement. Non-linearities in tip displacement are compensated 
by varying bourdon geometry or by adjusting the movement. 

4. Temperature Effects 

One of the most difficult problems that manufacturers of precision pressure 


gauges must overcome is the effects of temperature on the elastic element. Tip 














displacement is very susceptible to changes in ambient temperature. The temperature 
characteristics of a Bourdon depend on geometry and material. Fluctuations in 
temperature affect the elastic modulus, or stiffness, of the bourdon. Higher temperatures 
correspond to lower stiffnesses which mean greater tip travel. The opposite is also true. 

Thermal length changes affect the entire measuring system. The easiest and 
cheapest way to compensate for temperature effects is to manufacture elastic elements 
from materials having low coefficients of thermal expansion. The material most 
commonly used for this purpose is Ni-Span C. Although less affected by temperature, 
Ni-Span C does not completely eliminate the effects of temperature. Manufacturers also 
compensate for temperature effects by making all of the components of the measuring 
device from the same material. This approach eliminates any tip deflection relative to 
the movement caused by temperature changes. [Ref. 2,4] 

Some gauge manufacturers maintain gauge accuracy across fluctuating 
jeniperatiies with the addition of a temperature compensator within the movement 
mechanism. The device reacts to the same temperature that affects the bourdon. As the 
temperature changes, the device adjusts a setting within the movement to offset the 
temperature effect. Each compensator is custom manufactured to match and offset the 
temperature characteristics of the elastic element with which it is used. [Ref. 2] 

5. Technical Publications 

Only one publication was found that presented an analytical solution for bourdon 
behavior. In 1946, Alfred Wolf wrote An Elementary Theory of the Bourdon Gage which 


presented theory based on two elements of strain: strain due to “the bending of the walls 














in a transverse section through the tubing” and to “longitudinal extension and contraction 
parallel to the axis of the tubing”. In this publication, Wolf derived formulas for what he 
called the sensitivity (6R/R) and the torque (7) of the bourdon where dR was defined as 
the change in radius of curvature and torque was defined as the external bending moment 
required to restore tip deflection to zero for any given pressure. Based upon his 
theoretical results, Wolf concluded that the theory of pure bending in thin shells did not 
adequately model the behavior of the bourdon for the case a‘/R’P much greater than 
unity and that gauge sensitivity was inversely proportional to the first power of ¢ for thin- 
walled tubes. Wolf's tabulated results for “observed” and “calculated” sensitivities and 
torques for a steel elastic element were quite comparable to on another. 

6. Manufacturers’ Literature 

The logical progression in the search for literature was contacting companies that 
manufacture the elastic elements addressed and analyzed in this paper. A representative 
from the Stewart Warner Instrument Corporation commented that their gauges cover a 
very broad range of pressures and are “old standards” based on tried and true designs. 
On the rare occasion that a customer requests a pressure gauge not previously 
manufactured, the company has a “formula” for estimating the new design based on 
established gauge designs. The company representative also said he had no knowledge 
of an analytical model to determine tip travel. 

Excerpts from a pressure measurement handbook produced by the WIKA 
Instrument Corporation provided a brief history, the principles of operation, some design 


considerations, the effects of temperature, and load constraints of bourdon gauges. The 
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excerpts also presented equations for tip travel, relative change in curvature, and 


maximum bending stresses of tubular elements but, for proprietary reasons, all were 
expressed as functions of undefined variables. 

The Heise Division of Dresser Industries also provided some literature on dial 
gauges. Their biggest contribution to this literature survey was information concerning 
temperature compensation, methods of which were previously discussed. Like the 
Stewart Warner Corporation, Heise’s gauges are well established and rooted in 50-year 
old designs. Heise is, however, currently developing its own analytical model for 
Bourdon elements. 

The efforts of a literature survey yielded relatively little information concerning 
analytical studies of Bourdon pressure gauges. Very few technical papers have been 
published which analytically describe Bourdon theory. As a result, gauge manufacturers 


typically rely on experimentation and established designs when creating new gauges. 


C. RESEARCH OBJECTIVES 

The objectives of this research were to develop an analytical solution for tip 
travel, specifically in C-type Bourdons, and to compare results from the analytical model 
to results from finite element models and data collected from laboratory experiments. 
Because pressure gauge manufacturers most commonly rely on experimentation for the 
design of new gauges, a realization of the stated objectives could prove to be quite useful 


to them or to anyone interested in the applications of Bourdon theory. 
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Il. ANALYTICAL ANALYSIS OF A BOURDON TUBE 


A. NOMENCLATURE AND ASSUMPTIONS 
1. Nomenclature 
For clarity, symbols related to Bourdon geometry are annotated in Figure 4 and 
defined as follows: 
a - semi-major axis of the elliptical cross-section 
b - semi-minor axis of the elliptical cross-section 
t - wall thickness of the Bourdon element 
R - radius of curvature of the Bourdon element 
“- angle describing the length of the Bourdon tube 
/ - circumference of the Bourdon tube 
p - radius of curvature of the cross-section of the Bourdon tube 
¢ - angle describing the position of a point on the circumference of the 
Bourdon element 
k - complete modulus of elliptic functions 
K - complete elliptic integral of the first kind 
E - complete elliptic integral of the second kind 
E" - Young’s modulus 


E; - elastic modulus from the theory of thin shells 














Figure 4. Bourdon length segment and cross-sectional view 
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2. Assumptions 


In the process of developing this analytical model for Bourdon behavior, the 
following general assumptions were made: 
a) the circumference / remains constant during cross-sectional deformation, 
b) the length of the tube, defined as 4R, remains constant as the elastic element 
straightens, 
c) the wall thickness ¢ remains constant as the cross-section changes, 
d) the change in the angle ¢ (4 - ¢) is so small for small deformations that ¢ 
essentially remains constant during deformation, and 


e) the material of the elastic element is homogeneous and isotropic. 


B. MATHEMATICAL BACKGROUND 

As a preface to this section, detailed derivations of all equations are presented in 
Appendices A-D. For the sake of clarity and brevity, only simplified expressions are 
presented in the body of this paper. 

Some background concerning Jacobian elliptic functions and integrals is 
necessary for following the mathematical progression in this thesis. The following 


definitions 














and 


vl a 


E= j1—K-sith-6 do 22 


are complete elliptic integrals of the first and second kind, respectively, [Ref. 9, 10] and 
appear frequently throughout the problem formulation. The variable k , defined by the 


expression 
b 2.3 


on the interval 0<k <1, is known as the modulus of the Jacobian elliptic functions. 
Identities containing the elliptic functions sn(u), cn(u), and dn(u) are presented in 
equations 2.4.a - 2.4.f which reduce to the trigonometric functions sin(u), cos(u), and 1 


when k is zero. Elliptic function identities relevant to this problem formulation are: 


sn’ fu) + cn*(u)=1 2.4.a 
dn*.(u) + k* sn? (u)=1 2.4. 
sn(0)}=0 2.4.¢ 
en(0)=1 2.4.4 
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dn(0)=1 24. 


en(K)=0 2.4.f 


which are also defined in Refs. 9 and 10. Solutions for the following specific integrals 
are derived here for later use in the computation of strain energy in elastic Bourdon 


elements due to cross-sectional deformation: 


























s E(2-k* 
| Pa Ne ee ZEQ-K) y 2.5 
0 s(1-k)| a-k) 
K 2 2 
| ac *(u) do ices C2, Jel 32 i 
“0 wa-kKy | a-k) 5(1- ky 
= 24-(2 - ky" 5 2-E-(2-k*) 18-E-(2 - k*) 
dn ®.(u) du= a - | | - | - 
“0 15(-k*Y ad -k’) aa-k) 35-(1 - k’) 
27 
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C. PROBLEM FORMULATION 

1. Geometric Relationships 

Assuming that the volume of the Bourdon changes as deformation occurs, the 
change in the semi-major axis da is expressed as a function of the change in the semi- 


minor axis db by the equation 


-b-(E - 
wtb(E-K) 
Ea’ -K-b? 


da 2.8 


The derivation of equation 2.8 is presented in Appendix A. For a unit length, the volume 


of the Bourdon is defined by the equation 


Ven-ab 2.9 
The change in volume then becomes 
dV=x (a db + b-da) 2.10 


or 


(a? + b?)-E - 2Kb? 
Ea’ - Kb” 


dV=na -db 2.11 











by making the substitution for da in equation 2.10. 
Bourdon tip travel is a function of the sectional deformation that occurs as a 


result of pressurization. For an applied internal pressure, a curved circular tube will not 
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straighten. But a Bourdon straightens as its elliptical cross-section becomes more 
circular. Because of this behavior, it is important to understand how the radius of 
curvature of the cross-section, p, changes. As can be seen in Figure 4, the radius of 
curvature for a circular section is constant but changes as a function of the angle ¢ for an 
elliptical section. A simplifying assumption here is that, as a point is mapped from the 
elliptical shape to the more circular shape, the defining angle ¢ is approximately the 
same before (¢,) and after (¢,) deformation for small tip travel and minimal sectional 
deformation. 

Whereas p defines the radius of curvature, //p defines the curvature of the cross- 


section. The simplified expression for curvature may be written as 


b 
2:12 


1 
P | a2-(1- ksit-o)” 





The change in curvature, determined by taking the derivative of //p only with respect to 


dimensions a and 5, is defined by the expression 


(-a' - b’)-E 
a age 


K }-(1 - k’ sin’ ¢) +3-E| db 














a(4 ee ee! 
PL Bk sin?) ”- (Eat - Kb*) 

213 
The derivative of //p was not taken with respect to the angle ¢ because, as stated earlier 
in the assumptions, ¢@ is assumed to remain constant throughout pressurization and 


subsequent deformation. 














2. Incremental Strain Energy Due to Cross-Sectional Deformation 
Having derived the fundamental geometric relationships necessary to develop the 
analytical model, energy methods were applied to obtain a solution to the problem. In 


general, the elastic strain energy U in a beam due to bending is written as 


Us —— dx 2.14 


MsE I a() 2.15 


Figure 4 shows that, for an incremental slice of the cross-section, the wall of the tube can 
be modeled as a curved beam where d/(//p) is the change in curvature of the cross- 
section. From the theory of thin shells, the effective elastic modulus £, is defined by the 


expression 





E= ; 2.16 


where E’ is the elastic modulus of a straight beam. The area moment of inertia / is 


defined for a unit length slice of the tube by the expression 


_! 2.17 
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Rewriting the strain energy equation as symmetric in four quadrants, the 


expression for the total incremental strain energy due to cross-sectional deformation is 


Substituting equations 2.15, 2.16, and A.14 into equation 2.18, the incremental strain 


energy becomes 





Tt 
2, oe 
dU = a a: a") 1-Kk-sitt- do 2.19 


Pp 


and eventually, after considerable simplification, 


K 





2-a-E'E-C,-db" 
dU = 


x — Cydn*(u) + Cyd ®(u) + Cedn ®(u) du 
i-v) 0 





2.20 


where the constants C, through C, are defined as follows: 


4 
Cee se 
a’ (Ea’ - Kb’) 




















oath : 
C= E-K 2.20.b 
b* 
2 2 
-a-b 
cae | : -E- 3 2.20.c 
b 
C =9-E? 2.20.d 


The variables K and EF represent the complete elliptic integrals of the first and second 
kind, respectively, and are defined in equations 2.1 and 2.2. The integrals of dx? (u), 
dn® (u), and dn* (u) are defined by equations 2.5, 2.6, and 2.7, respectively. 
3 Incremental Strain Energy Due to Bending 

To determine Bourdon incremental strain energy due to bending, or straightening 
in this case, an expression was developed for the total bending moment due to an applied 
internal pressure. Referring to Figure 5, the net moment M is found by summing the 
moments acting on the end cap and walls of the tube about point O. The details of this 


derivation are shown in Appendix C and the resulting expression is 


—_ 


=-P-aR-D-(1 - cos@)-k 221 


where P is the applied pressure, R is the radius of curvature measured to the centerline of 


the tube, and D is a constant for specified dimensions a and 5 defined as 
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end cap; Bourdon tip 


Figure 5. Bourdon length segment and cross-sectional view depicting the elemental area 


for the derivation of the net moment acting on the elastic element 
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Again, K and E represent complete elliptic integrals of the first and second kind, 
respectively, as defined in equations 2.1 and 2.2. Substituting equations 2.21 into 
equation 2.14, the total strain energy due to bending of a length of tube that subtends the 


angle with the center of curvature becomes 


(-P-a-R-D-(1 - cos(B))) p 48 2.23 
2E'I 





“0 


where the moment of inertia J is now defined as 


Ie=-((a+ t)(b +t) - ab*) 2.24 


The incremental strain energy due to bending was found by differentiating U, with 


respect to a, b, P, and R as follows: 


d d d d 
dU, =— Ur fine UL prigpeees U, + U, 2.25 


da db dP dR 


To simplify the results of this lengthy equation, the incremental strain energy due to 


bending is rewritten as 
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dU,=A-(B +Z) 2.26 


where A, B, and Z are defined as follows: 





P-a’-R?.-D 
As 6-2 - 8 sinf a) + sin(2: 
ee (6-4 (a) (2-2) ve 
pee RP E-D id ppl aeprRri dpa noe 
(Ea’-Kb") \ da db an 
and 
Z=(R-P-dP) + (1.5-P-D-dR) 2.26.¢ 


According to the principle of virtual work defined in Ref. 11, the change in 
external work is equal to the change in elastic strain energy (dW = dU) if a configuration 
satisfies geometric constraints and remains continuous. Assuming that the Bourdon 


satisfies these conditions, 


dWedU_ + dU 2.27 
xs b 


where the incremental external work acting on the tube is defined as 


dWe=P- dV 2.28 


or 
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(P+ dP)-dVedU, + dU, 229 


Making the appropriate substitutions into equation 2.29, the resulting equation is written 


as 


(a +b)-E-2K-b 
Ea’ -K-b" 


(P+ dP)-x-a db= 





2aEIC,-db" K 
—________—- (c -dvtuse Cdn ue dau) du + 
2 2 3 4 
l-? 9g 
P-a-R-D 





(1- v’)-(6p — Sins + sin(2-p)) x 
4E'-] 





PREE-©) 4 ppl aepr(o pl abePRaper a 
(Ea’-K-b4) | da db : 





2.30 
which is quadratic in db where dR is also unknown. See Appendix B for derivation 
details. 

4. Tip Travel 
The equation for tip travel was derived by making the assumption that the length 
of the Bourdon remains constant for any given pressure and by applying Castigliano’s 


second theorem. The assumption for constant length, L, written as 
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wRaL 231 


implies that 
uw dR + R-du=0 252 


or that incremental tip travel dR can be calculated using the expression 


dpe 2.33 


ul 
Whereas R and yz are known quantities, the variable dy is undefined. Referring to 
Figure 6, 


duz=- 6G; 2.34 


by geometry where @,, the angle of rotation of the tip due to the applied dummy moment, 


is described by the equation 





which is the application of Castigliano’s second theorem [Ref.11]. When the dummy 


moment, M74, is applied, the equation for the net moment becomes 


3 
M=-P-aR-D-(1 - cos§)-k ~ M,k 2.36 














Hi > Hy 








uR = Constant 
dR=R,-R, 


[es =[s- 


Figure 6. Bourdon length segment and tip deflection for the application of 


Castigliano’s second theorem 
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Substituting equation 2.36 into equation 2.35 and setting the dummy moment equal to 


zero, the expression for tip rotation is obtained by integrating 





In this formulation, tip displacement is calculated iteratively. For this model, knowing @ 
defines du which then allows the determination of incremental tip travel dR. Once @R is 
known, db is found by solving equation 2.29 and da is calculated using equation 2.8. 


Subsequent values of dR are then calculated using incremented values of a and b. 


D. SIMULATION PROGRAM AND ANALYTICAL RESULTS 

1. Source Code Description and Limitations 

Considerable effort was expended in expressing all of the equations in Chapter II 
as functions of a, b, k, EF, and K to simplify the programming portion of this research 
effort. Having derived the requisite equations to determine tip travel and sectional 
deformation, a program was written in FORTRAN to iteratively calculate solutions for 
user-defined Bourdon properties, geometries, and operating pressures. The source code, 
listed in Appendix E, defines all program variables, prompts the user for input, and then 
calls seven subroutines to iteratively calculate tip travel and corresponding changes in tip 
geometry. 

The topics of temperature effects and the standard methods of compensation 


discussed in Chapter I were not addressed in this attempt to develop an analytical 
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solution for Bourdon tip displacement. The analytical model presented in this paper only 
calculates tip travel and cross-sectional deformation as functions of pressure. 

2. Analytical Results 

Execution of the source code Program Bourdon was very sensitive to the 
geometric relationships between a, b, and R and worked best for a:b ratios near 2:1 and 
for 0 < R< 34 mm (0 < R< 1.34 in). The program ceased to work beyond one iteration 
once R exceeded 33 mm. 

The following table indicates analytical results for fictitious Bourdon elements 
each fabricated from phosphor bronze having an elastic modulus of 110 GPa (16 x 10° 


psi) and a Poisson’s ratio of 0.34: 


Table 1. Analytical Solutions Generated from Program Bourdon 


Bourdon Radius of Wall Angle Max. Tip Travel 
Element Curvature | Thickness | Describing | Applied (dR) 
{R) (t) [mm/in] Tube Pressure [mm/in] 
{mm/in] Length (1) (P) 
[deg/rad] | [MPa/psi] 


10.000 / 30.000 / 0.300 / 250 / 1.00 / 
0.393 1.181 0.012 4.363 145.00 


10.000 / 30.000 / 0.100 / 250 / 1.00 / 
0.393 1.181 0.004 4.363 145.00 
10.0 / 30.000 / 0.300 / 330 / 1.00 / 


0.393 1.181 0.012 5.760 145.00 


0.300 / 250/ 
0.012 4,363 
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14.000 / 30.000 / 0.100 / 250 / 1.00 / 112.41 / 
0.551 0.197 1.181 0.004 4.363 145.00 4.43 


5 10.0 / 50.000 / 0.300 / 7.12/ 
0.393 1.97 0.012 4.363 0.28 

7 14,000 / 0.100 / 250 / 1.00 / 33.31 / 
0.551 0.004 4.363 145.00 1.31 






20.000 / 


0.787 






7.500 / 0.100 / 250 / 1.00 / 4,32 / 
7.500 / 0.100 / 250 / 1.00 / 14.57 / 
1.181 0.004 4.363 145.00 


Looking at the results of the analytical model, some of the tip deflections are 





excessively large while others are of more reasonable magnitudes. These extremes in 
output are currently unexplained and under investigation. While the validity of the 
analytical model’s results is unknown, an important observation is made concerning the 
trends of the output. As stated previously in Chapter I, Bourdon elements having flatter 
cross-sections and thinner walls exhibit larger tip displacements than rounder, thicker 
elements for the same applied pressure. The results presented in Table 1 show trends 
that agree with Bourdon displacement theory. 

Comparing Bourdon elements 1 and 2 from Table 1, element 2 has a thinner wall 
than sleitient 1 and consequently experiences significantly more tip travel. Looking at 
the major axes of elements 2, 6, and 9, element 6 has the flattest cross-section followed, 


in increasing order, by elements 2 and 9, respectively. As expected, element 6 exhibits 














the greatest tip movement and element 9 exhibits the least when all other parameters are 
held constant between the three specified elements. Similar comparisons can also be 


made between other elements in Table 1. 
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Ill. EXPERIMENTAL AND NUMERICAL ANALYSES OF A 
BOURDON TUBE 


A. EXPERIMENTAL ANALYSIS 

Because analytical studies are based on assumptions, it is good practice to 
compare analytical results with experimental data. Doing so verifies the validity of the 
mathematical model and its underlying assumptions. For the model developed in this 
thesis, a good correlation between analytical results and experimental data implies the 
successful application of this analytical model to future designs of Bourdon gauges. 

The experimental results presented in Tables 2 and 3 were obtained by two naval 
facilities: the calibration laboratory at the Regional Repair Center (RRC), Norfolk, and 
the Mechanical Instrument Repair Calibration Shop (MIRCS) in USS SHENANDOAH 
(AD-44), homeported in Norfolk. The requested information was faxed to SIMA, 
Norfolk and was then disseminated to the aforementioned facilities. Pertinent 
dimensions, elastic element information, and Bourdon tip travel results from the 


experimental studies are: 














Table 2. Experimental Data from RRC 
























Radius of 





Tip Travel 
(dR) 
[in/mm] 


Gauge Bourdon Applied 









Pressure (P) 
[psi/MPa] 


Curvature 
(R) [in/mm] 





Company Material 





and Range 
[psi] 


Weksler Alloy Steel 0.5737 / 0.2220 / 1.361 / 300.00 / 0.25 / 
ne | em | "om | sane | oe |e 
Mensor Stainless 0.5559 / 0.2936 / 3.0843 / 1000.00 / 0.4375 / 
Heise Stainless 0.4108 / 0.2758 / 3.5129 / 3000.00 / 0.25 / 


Table 3. Experimental Data from MIRCS 




























Radius of Tip Travel 
Curvature (R ) (dR) 
[in/mm] [psi/MPa] [in/mm] 


Stainless 0.5 / 0.115 / 0.771 / 500 / 0.131 / 
Beryllium 0.5 / 0.190 / 1.033 / 600 / 0.129 / 


One key parameter missing from Tables 2 and 3 is the Bourdon wall thickness of 


Applied 


Pressure (P) 






Major Axis (2a) 
{in/mm] 


Minor Axis (2b) 
{in/mm] 


Bourdon 
Material 















each gauge. Whereas the major and minor diameters were easily measured, wall 
thicknesses were neither measurable nor available from the respective manufacturers due 


to the proprietary nature of the information. 
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B. FINITE ELEMENT ANALYSIS 


One rapidly growing approach to engineering problem solving is the numerical 
technique of finite element analysis. These techniques (1) discretize a domain with a 
finite number of finite elements and a finite number of nodes per element, (2) introduce 
nodal points and appropriate shape functions to represent the trial solution in each 
element, and (3) transform the resulting set of differential equations into a series of 
matrix equations to change a system from an infinite degree of freedom system to a finite 
degree of freedom system of appropriate size. 

Several powerful finite element software packages exist today that compute 
solutions for complicated engineering problems. For the purpose of this research, the 
ability to build, constrain, and pressurize a Bourdon tube using computer graphics greatly 
enhanced the study of Bourdon behavior. As powerful as this tool can be, however, this 
type of analysis did not preclude the effort to conduct laboratory experiments but rather 
supplemented the study of Bourdon behavior. Once the models were constructed and 
drawn to match the dimensions and material properties of the Bourdons used in the 
experimental studies, tip displacement results were quickly determined. 

The finite element solutions presented in this paper were generated using I-DEAS 
Simulation software, software that builds geometry and finite element models, analyzes 
the models, and evaluate the results. Each of the finite element models presented in this 
paper was constructed using the following parameters: (1) major axis 2a, (2) minor axis 
2b, (3) wall thickness ¢, (4) radius of curvature R, (5) Bourdon arc length y, (6) elastic 


modulus E, (7) Poisson’s ratio v, and (8) applied pressure P. 


3D 











Assumptions for wall thickness and arc length were made based on design 
information presented in Refs. 1 and 5. Boundary conditions were applied at one end of 
each Bourdon model. The socket or pressure end was clamped to prevent rotation and 
translation in the x, y, and z directions. The tip end of each model was sealed but 
otherwise left unconstrained to allow maximum freedom of movement. Internal 
pressures corresponding to experimental test pressures were applied to the walls and end 
cap (tip) of each finite element model. 

The numerical models were analyzed using linear, thin-shelled elements whose 
thicknesses were modified to match the estimated wall thicknesses of the Bourdon 
elements listed in Table 1. The quality of the mesh for each element was established 
using criteria for distortion. On a scale of -1.0 to 1.0, an element is perfectly square for a 
distortion value of 1.0. Meshes exhibiting distortion values less than zero are very bad 
and should be avoided whereas distortion values greater than 0.4 are acceptable. In 
modeling the Bourdon elements using I-DEAS, all of the meshes exhibited distortion 
levels greater than 0.4. The elements that were distorted, however, consistently occurred 
along the length of the element at the major axis. 

The following table summarizes dustherieal solutions for tip displacement in 
Bourdon elements having the same geometries and material properties as those presented 


in ChapterI1.Table 1: 
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Table 4. Numerical Solutions for Table 1 Bourdons 

























Radius of Wall 





Bourdon Angle Max. Tip Travel 


























Element Curvature | Thickness | Describing | Applied (dR) 
(R) (t) [mm/in] Tube Pressure {mm/in] 
{mmvin} Length (11) (P) 









[deg/rad] | [MPa/psi] 
10.000/ 5.000 / 30.000 / 0.300 / 250 / 1.00 / 3.26 / 
ae ee 
14.000 / 5.000 / 20.000 / 0.100 / 250 / 1.00 / 42.8 / 
pM eo ee | | tee | 
7.500 / 5.000 / 20.000 / 0.100 / 250 / 1.00 / 8.50 / 
LT [ee [toe | tr | te | so | sae | 





Bourdon elements 1, 7, and 8 were chosen from Table 1 at random for modeling 
in ]-DEAS. As for the case of the analytical model, the validity of the magnitudes of tip 
deflection generated by the finite element technique is also uncertain without 
comparative data from gauge manufacturers or more extensive experimental work. 

Comparing the results from Tables 1 and 4, the numerical solutions follow the 
same trends as those from the analytical model but at greater magnitudes of tip 
deflection. Tip travel is 53%, 22%, and 49% greater for I-DEAS elements 1, 7, and 8, 


respectively, than for the same analytical elements in Table 1. 
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IV. CONCLUSIONS AND RECOMMENDATIONS 


A. CONCLUSIONS 

The objective to analytically derive a model to predict Bourdon behavior was a 
very challenging one. The very mechanism by which Bourdons move, as described in 
Chapter I, indicates the level of complexity of this engineering problem. The simple 
beam and thin-shell theories and equations used in the development of this analytical 
model for tip travel in a Bourdon tube have produced solutions whose validity is yet to be 
determined. 

The task of comparing analytical, numerical, and ees results was 
hampered by insufficient experimental data, questionable numerical solutions, and forced 
assumptions concerning Bourdon geometry and material. The fact that Program 
Bourdon produced output that followed expected trends is encouraging. The fact that the 
code runs for more than one iteration only for certain Bourdon geometries warrants 
continued investigation. 

The experimental data collected from RRC and MIRCS was very limited due, in 
some part, to the underspecification of the requested information. Ideally, the 
dimensions of the elastic elements should have come from the gauge manufacturers, but 
as stated previously, that information is closely guarded by the companies that design and 
manufacture the gauges. In addition, the tip travel for each bourdon element should have 
been measured over the operating range of the gauge and not just at the maximum 


allowable pressure. 











The accuracy of the numerical results generated in I-DEAS was questionable due 
to the excessive magnitudes of tip deflection for the applied pressures. Because the 
experimental and numerical results for tip travel were each somewhat inconclusive and 
because the analytical results were constrained by the geometric relationships of 
dimensions a, b, and R, the comparative study of results yielded little useful information 
for us to be able to state with conviction the accuracy of the analytical results. 


Consequently, the validity of the analytical model remains under investigation. 


B. RECOMMENDATIONS FOR FUTURE RESEARCH 

The research results presented in this paper are rather limited implying that a 
considerable amount of work still needs to be done to develop an acceptable analytical 
model for Bourdon behavior. The following suggestions are made to continue the 
research effort: 

1. Peruse Chapter VII (Tubular Manometric Springs) of Ref. 12 for its detailed 
discussion of Bourdon tubes and compare the assumptions and derivations presented in 
this thesis to the theories presented in Elastic Elements of Instruments. This book is 
highly recommended by the Heise Division of Dresser Instruments as an excellent 
resource in developing an analytical model for Bourdon behavior. 

2. As a means of eliminating the geometric constraints of the analytical model, 
modify the model using less restrictive geometric assumptions by independently 


considering variable Bourdon length, circumference, and angle ¢. 
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3. Collect more thorough experimental data covering a much broader operational 
range of pressures to better analyze and predict Bourdon behavior and eliminate as many 
assumptions as possible concerning Bourdon geometry by obtaining comprehensive 
measurements. 

4. Enhance the analytical model to include a stress analysis to predict failure 
from fatigue due to cyclic loading or from overpressurization and exploit the ability of 
finite element software packages to calculate stresses and strains within finite element 


models for comparison purposes. 
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APPENDIX A. GEOMETRIC RELATIONSHIPS 


Geometric relationships are derived in great detail in this appendix. These 
relationships ultimately serve as the building blocks for the formulation of the analytical 
model for Bourdon behavior. Derivations for da, the change in the semi-major axis, dV, 
the volumetric change of the tube of unit length, and d/1/p), the change in the curvature 
of the cross-section, are presented in this appendix. 


Referring to Figure 4, the equation for the elliptical cross-section is 
2 
Zz —— 
a 
and the area, A, is 
A2nt a: b A.2 
The volume of the elliptical tube with length unity is 
Venrab A3 
and the change in volume is 


dV=n-(a db + b-da) A4 


The circumference, 1, of an ellipse is defined as 














1=4. ds A.5 


where ds is an incremental length of the circumference defined by the radius of curvature 


of the cross-section p and the angle ¢@ Referencing Figure 4 for small d¢, 


ds= fd? +z A6 


For 


x=b-cos¢ A7 
and 
za: sing A8 


the following expressions are true: 


dx’=b* sin’ ¢- dg” Ad 


and 


dz*= a? cos’. $- dg” A.10 


Substituting equations A.9 and A.10 into equation A.6 and dividing both sides by the 
square root of a’, 
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b 
ds=a- — sin’ + cos’ ¢-dg All 
a 


Substituting in the trigonometric identity cos’¢ = | - sin’ ¢ and rearranging terms, 





2 


Masha yo A12 
a 


rn, 
Lom 


apse A.13 


on the interval 0< & <1. The incremental arc length ds is now written as 


ere, ~ k’ sin’ ¢- dg A.14 


so that the circumference becomes 


A 
2 


I=4-a 1-k’* sin’ $ dé ATS 
“0 


and then 


l=4-a-E 

















where E is the complete elliptic integral of the second kind defined by equation 2.2. For 
the assumption that / is constant, d/ = 0. Taking the derivative of equation A.16 and 


equating it to zero yields the following expression: 


rake 


Solving for dk, the derivative of both sides of equation A.14 is taken and written as 


+ E-da=o A.17 








k-dk=—?.(a db ~ b-da) A.18 
a 
then as 
dk= p -(b-da - a db) 
kw A.19 


and eventually as 





b*.da- ab-db 
dk= A.20 
a’ Ja - b* 
after making the substitution for k defined by equation A.13. From Ref. 9, 
d E-K 
pais A21 
dic k 
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where K is the complete elliptic integral of the first kind defined by equation 2.1. 


Substituting equations A.20 and A.21 into equation A.17 yields 
E-K 
b | } 
k 


which simplifies to 


b-da - a-db 


ka 





+ E-da=0 A.22 








(E-b? - K-b")-da + (K-ab - Eab)-db + (Ea? - Eb!)-da_ 


A.23 
a’ - b’ 
and eventually to 
Ea - K-b’ 


the ultimate expression of the change in the semi-major axis written as a function of the 
change in the semi-minor axis. Substituting this expression into equation A.4, the 


volumetric change in the Bourdon may be written as a function of db as follows: 














den ehh a A25 
Ea -K-b" 
After some algebraic manipulation, the change in volume is written as 
2 2 3 
dVenr- Peat eels db A.26 


Ea’ -K-b? 

















Referencing Figure 4, the radius of curvature p of the cross-section is defined by the 


relation 


A.27 
dg 


1 
Fz) dz 2] A.28 
1+ {— 
dx 
The objective now is to write curvature as a function of the parameters a, b, and ¢. 


Recalling equation A.1, its derivative is 








xdz zdz 
pea + eras! A.29 
a 
Solving for dz/dx, 
2 
ae -ax A30 
dx pf, 
so that 
dz\? ats? 
& “Ete fe 


The next objective is to write dz’/dx’ in terms of a, b, and ¢. By definition, 
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dz d fdz a 
dv? dx de 





dz do ~al-x 

de de b*z A.33 
or 

diz _-a do x 

“a ee A34 


which, when the derivative is taken, becomes 


dz 

~] = aia 

diz -a’ oe ( - 
——a—.} A.35 


dx? ? Zz 


Substituting equation A.30 into equation A.35 and simplifying, equation A.35 may be 


rewritten as 








and eventually as 














4 
-a 


23 
bz 


A.37 





d. Z_ 
ae 
Now substituting equations A.30 and A.37 into equation A.28, curvature is expressed as 


Pee ee eet: A.38 
p aa 2)" 


a ae | A.39 


Substituting equations A.9 and A.10 into equation A.39, curvature is written as 


1 -a'-b* 
Pz 4 2.3 4,2 2 ,,15 A.40 
P (b’-a’-sin'-¢ + a-b’-cos’$) 
OT as 
1 -ab 
A.41 


© (b* sin? 6 + a-cos’-)” 


Recalling the trigonometric identity cos’¢ = / - sin’¢, and dividing the numerator and 


denominator by a’, 
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1 -b 
> |e 1s A.42 
— sin’ ¢ +(1- sin” $) “a 
Z 
and ultimately 
Lis: -b 
PO hc - k? sin?) A.43 
Considering the absolute value of the curvature 
1 b 
a. A.44 


e Ba kesin? gy” 


which now correlates to the definition of curvature for an elliptical cross-section as 


defined in Ref. 13. The change in curvature may now be written as 
d- () =—-do A.45 


where 


fe. see 
o af(1—k sin’ gy” 


A.46 


and 











do= fas 5a a4 -dk A.47 
da db dk 7 


The derivative is taken with respect to & because k is a function of the dimensions a and 5 


as defined by equation A.13. The partial derivatives of p with respect to a, b, and k are 


2: : 

f p\ daw (ia 4 ging) da oor 
da b . . 
2 a\t6 a (1 - ke sin?) db 

—_— . =-—: ~ *$ . 5 

db 4? A.47.b 

2 - 3 
eae a oer Ate 


so that equation A.47 becomes 


‘ 2 -ae-Je- dle: sin’: 
doz ra i <a} = ik? sin? 6)” “2 os I ~ k?. sin’ ¢ 


b b? 
A.48 
Recalling equations A.13 and A.19, the expression kdk may be written as 
k-dk= i b-da - a db 
ae (b-da - a db) A.49 
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Making the substitutions for da and kdk and simplifying, equation A.48 becomes 


(2b? — a“)-E - Kb? 


(1 - ke sin? g)'”- db - 
Ea’ - Kb" 


af 
acer 








ea 
b* - a 


3-E-sin® ¢- 
Ea’ - Kb’ 


A ~ k* sin’. ¢- db 





A.50 


which is even further simplified by manipulating the following terms: 


sin’. $-aft - Pst gm ~ 1 sin? g-((1 - sin? 6) - 1) A.51 


k 


or 


sin’. ¢-,/1 - sat ge | ~ sin? gy” — ft - 2. sn? | A.52 
k 


Making the substitution of equations A.13 and A.52 into equation A.50 , the expression 


for the change in the radius of curvature becomes 


2 
Poesy | ~ k* sin’. ¢- db + 


Ea’ - Kb’ 


(1-2 sin? g) > db 


a (2b'-a)-E-Kb  3a'E 
b’ Ea’ - Kb* Ea’ - Kb" 


A.53 


which, with a little more manipulation, becomes 


a3 














3 2 2 2 2.2 ,.15 
geese ie ype EE Re ee 
Ea’ - Kp? b? 


A.54 
Now substituting equations A.46 and A.54 into equation A.45, the change in curvature is 


written as 


(a + b*)-E 
seers r 


3B 





K 











2 
a(t) a ie i sin?) s 
PL Bc. - ke sin’ $)"”- (Ea? - Kb?) 


A.55 
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APPENDIX B. INCREMENTAL STRAIN ENERGY DUE TO 
CROSS-SECTIONAL DEFORMATION 


In general, elastic strain energy for a beam subjected to transverse loading is 





expressed as 
y 2 
Us uM dx B.1 
2-ET 
“0 


So as not confuse the symbols for elastic modulus and elliptic integrals of the second 
kind, elastic modulus for a straight beam will be designated E’ , elastic modulus derived 
from the theory of thin shells will be designated E, , and elliptic integrals of the second 
kind will be designated F for all remaining calculations. For an initially curved beam, 


the moment Mis defined by the expression 


1 
MsE -I- B.2 


M B3 





-0 


where / is now defined as 
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MsE I a(t) B.4 





B. 2 
gq wet (2) . BS 
xs 2 


I=— B.6 


and ¢ is the wall thickness of the Bourdon element. The integral term is found by 


squaring equation A.45 which results in the expression 


2 


-(1 - k* sin? 6) |-db 

















-b* 
2 22 425 po 2 2 . 
a-(i-k’-sin-¢) -(Ea - Kb‘) 


B.7 


which upon squaring becomes 
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2 ty, 2 
(a Ee 


2a 
ee 
b 


b? b? 9: ‘ 
er rs Gena - 
a (Ea ~ Kp) (1 ~ k* sin? 4) (1 ~ k* sin? gy (1 ~ Kk’ sin? 4) 
B.8 
For simplicity, equation B.8 js rewritten as 


Cha ag See 


meee Or ta ee es 2.3 sf 
(1 ~ k* sin? 4) (1 ~ ke sin? gy (1 - k* sin? 4) 


Where the Constants C, 


B.9 
through C, are defined as follows: 
4 
CF 4 ~ 2,2 B.9.a 
a-(Ea*- Kp ) 
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Rewriting equation B.5 in terms of the constants C, through C, and substituting 
equation A.14 for ds, the equation for incremental strain energy due to a change in the 


cross-section becomes 





dU. = 
XS 
zr 
2 
2-E'I C C C 
‘C-db* AL a OEE AE SIM RY PASSER | fs sin?. dp 
1-9 (1- sin? g)” ke sin’. gy’ (1 - sin’? g) 


B.10 


With further simplification, incremental strain energy due to sectional deformation is 





written as 
dU = 
xs 
v 
2 
2E' La 2 C, c, Me 
raha cell 2.2 ra 2.2 Tas 2.2 4.5 dp 
1-7 (1 - k’- sin’. ¢) (1 - k’- sin’ 6) (1 - k’- sin’: $) 


B.11 
All of the terms in the denominator of the integral in equation B.11 may be written in 


terms of elliptic functions by letting 


sin: ($)=sn- (u) B.12 


and 
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cos-(¢)=cn (u) B.13 


so that 


d(sind)=cos¢-d¢=d(sn(u))=cn(u)-da(u)- du B.14 
Solving for do, 
dd=dn(u)- du B.15 
In addition, 


1 - k?-sin® (¢)=1 - k* sn? (u) B.16 


where 


1 — k* sn®(u)}=dn’-(u) B.17 


by definition in equation 2.4.b. Equation B.11 is further simplified by evaluating the 
integral in parts and by making substitutions with equations B.15 and B.17. The first 


term in the integral is written as 


ra 
2 K 
(1-2 sin? gy? dg | (d® (uy)? dn (u) da 


+ » 











K EK 
=| (da*(u)-dn(u)) a=] da *(u) do B.18 


“D - 


Applying equation 5.131.3 from Ref. 14 for m = -4, 





K 
dn *-(u) du= 
“0 
. K K 
é - k?. do *. (u)-sn(u)-en(u) - 2:(2 - k’)- | dn *. (uy du + | udu 
3-(k) ‘0 o 
B.19 
Using the elliptic function identities defined in Chapter I and the relationships 
(cty’=1 - 2 B.20 
and 
K 
| (dn. (u)) du=E B.21 
“0 
as defined in Refs 9 and 10, equation B.18 is rewritten in simpler form as 
= K 222k} = 
do *(u) da=s-———_ + gees dn * (a) du B.22 
Jo eG sie) -atiek tg 


Again, applying equation 5.131.3 from Ref. 14 to solve for the final integral where 
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m = -2, equation B.22 is further simplified to the expression 











K 2 
- 1 2E(2-k 
| dn *.(u) du=——_ a -K B.23 
“0 30-kK)] a-k) 
where, for future reference, 

K 
aci?.(u) dm —= B.24 
9 1- k’) 


Still in the process of simplifying the equation for the incremental strain energy in 
Bourdons due to changes in cross-section, the next integral to be evaluated is 
w 
2 K 
(1 - Ke sin’ g)°? dpe | (da®.(u)) **- da-(u) du 


“0 ~0 


K 


K 
= | (do? (u)-dn-(u)) du=| dn © (u) du B.25 
“9 “OD 


Applying equation 5.131.3 from Ref. 14 for m = -6, 














K 


da * (a) dus 
~D 


1 





K K 
~ k?. da’*.(u)-sn(u)-en(u) - 4:(2 - k). dn *(u) du +3: dn *.(u) du 
5-(k) “0 “0 


B.26 


Substituting equations B.20, B.23, and B.24 and applying the elliptic function identities 


defined in Chapter I, 
: iy 4(2-k) 12 E(@-k) 3-E 
dn” (a) arash ane ee K}- ae B.27 
“0 ia-kKy | G-k’) x01-K) 





And finally, applying the same procedure to the denominator of equation B.11 


corresponding to the constant C, , 


wt 
2 


K 
(1 - kK sin? gy"? de | (da? (uy *”. dn (u) du 


0 - 


K K 
7 (da'°.(u)-da-(u)) dus] dn®(u) du 
“0 -0 


B.28 
where, applying equation 5.131.3 from Ref. 14 for m =-8, 
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K 


da ® (u) du= 
0 


K K 
7 k*. do * (uy): sn(u)-en(u) - 6: (2 - k*)- dn’ °:(u) du + 5: | dn *-(u) du 
7-(k') “0 “0 


1 





B.29 


Substituting equations B.20, B.23, and B.27 and applying the elliptic function identities 





defined in Chapter I, 
. 24-(2 - k*y" 5 2E-(2 -k*) 18-E (2 - k*) 
sito ae) MOWED, |__2_—_|PEG-¥) _) 2G 
“0 m5(1-k7) 21-1 ~ ky (1-k) 35-1 -k*) 








B.30 
The equation for incremental strain energy due to changes in cross-section may now be 


written as 


2aE' LC, db" K 


dU : : 
(i-¥) Jp 


xs 





~4 -6 -8 
| c, aa -(u) + C,-da®.(u) + C-de’®-(u) | da B31 


where the integrals of dn“ (u), dn*(u), and di® (u) are defined by equations 2.5, 2,6, and 


2.7, respectively. 
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APPENDIX C. INCREMENTAL STRAIN ENERGY 
DUE TO BENDING 


As defined in Appendix B, the total strain energy due to bending in a beam is defined as 





CA 


where M represents the net moment acting on the bourdon as a result of pressurization 
and E’ is Young’s modulus. The detailed derivation of the net moment acting on the 
elastic element represents a very large segment of this appendix. 

Referencing Figure 5, the elemental area of the bourdon tube is defined by the 


expression 


dA=(R + pcos® )-dé-ds C.2 
Recalling equation A.14, equation C.2 may also be written as 

dA=a-(R + pcosd )-./1 - k’-sitt--d0-do C.3 
The force acting on the elemental area is defined as 

=} 

dF=P-dA C.4 


where 


3 
P=P. ((cosa- cosé)-i + (cose: sin@)-j + sine-k) C.5 











The angle @ describes the direction of the normal component of pressure acting on the 


wall of the tube as shown in Figure 5. Substituting equations C.3 and C.5 into equation 


C.4, the elemental force acting on the elemental area of the bourdon is 


_~4 
dF=P.((cosa-cos@)-1+ (cose sin8)-j + snaw-k)-| a(R + @- cos) afl ak sin’. 4- dé 4s | 


where the corresponding moment arm to the elemental area is defined as 


—_ 


=-R-(cos8-1+ sinf-j) + (R + o-cosd)-(cos6-i+ sin?-j) + o-sind-k 


Relating the angle a@ to the angle ¢ in Figure 6, 





cosa = tl 
1- kK sirf- 
and 
sim = =) 


aj1- k’-sirt-6 


As a result, equation C.6 may be expressed as 


— 
dF=P-a-(R + pcos® )- 
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(cos-cos8)-i + (cosd-siO)-j + ee -do- de 
a 


C.6 


C.7 


C8 


C9 


C.10 








The moment acting on the elemental area due to an elemental force is defined by the 


equation 


dMe=dr = dF C.1] 


Temporarily factoring out Pa(R+pcos¢)ddd@ from equation C.10 for simplicity, the 


vectorial components of the resulting cross product of equation C.11 are: 
costs p' sin - —-sinp-(R-sir8 - R-sirB + psp) | C.1l.a 
a 





cosd-cos0- p- sinb - D cos0 - R-cosB + p-cosd-cos@) 
a 





‘J C.11b 


and 


((R-cos8 - R-cosB + p-cosb-cosO)- cosdsirB — cosbcosO:(R:sir8 — R-sirB + p-cosp-sirB))-k 
C.ll.c 
For the sake of brevity, the cosine function will be abbreviated by “c” and the sine 
function will be abbreviated by “‘s” from this point to the end of Appendix C. Continued 
simplification begets the following equation for the moment acting on the bourdon tube: 


21 B 
M= Pars pet) poh 0 - PR p:cd )-sb-sO + Presb do 
a a 
“0 0 
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20 B 
fa Par: p-co )- 


“0 “0 


pco-sb-cO - PR + pc )-so-c8 + PRB dO do 
a a 





20 B 
+ P-a-(R+ p-co)-R-cb-(c8-sB - sO-cB)-k dO do 
“0 ae) 


C.12 
Factoring back in the expression Pa(R+ pcos¢)d¢d@ and evaluating the x-component of 


the equation C.12 first, 


2% B 
P-a- (R + p-co)-| p-co-sO - Epes + bsp -sd dO db =0 
a a 


C.13 
because sing is an odd function on the interval 0 to 2x. Similarly for the y-component of 
equation C.11, 

2% B 
-P-a: (R + p-cob)-| p-cb-cO - aR + p-co)-c8 + © RB -sb dO db=0 
a a 
“0 “0 


C.14 


Evaluating the z-component of equation C.11, 
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2n [Bp 
P-a-R: (R+ p:ch)-ch-(cO-sB - sO-cB) dO dp= 
/0 0 


4P-a-R 
3 





(1 - eB): = HORS ®) 
k 








C.15 


Therefore, the final expression for the moment acting on the walls of the tube is 


2 
> 4P-a 


M= = 
3 








(- 6, + GE K) C.16 


k? 
Again referring to Figure 5, the net force acting on the end cap of the bourdon element is 


eee 
Fopt Pr abj CAT 
and the corresponding moment arm is 
— 
dr=R-(1 - cB)-i- R-sB-j C.18 


Performing the cross product, 


—S OS 


Map tt x dF pt PRE ab-(1 = c8)-k C.19 
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Therefore, the net moment acting on the tube is the sum of the moments acting on the 


walls and end cap of the Bourdon as shown in next equation: 


























——_—_— 3 Ul 4a/K - E 
M,47M +M, =P aR - oe a +(2E-K)]-ab}k 
C.20 

For simplicity, the new variable D is defined as 

Dea 3 Ee ORER) C.21 
such that 

— 

Moz aR D-(1 - cB)-k C.22 


Substituting equation C.22 into equation C.1, the total bending strain energy is expressed 
as 


# 2 
7 (P-aR-D- (i - c8}) R 


a6 
ELI C.23 


“0 


The incremental strain energy is calculated by differentiating U, with respect to the 
dimensions a and 5, the applied pressure P, and the radius of curvature of the tube R such 


that 
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de uaa Sule. d_y).ap + |4_u}-aR C.24 
da \db dP dR 
where 
2. p3. 
Wg) ae ES es ene Sagi ta Dep ds C.25 
da 4EVI da 
2 2 ps 
SAG hm eee ian eeny (ool a ee 
db 4E'I db 
d P-a-R*-D? 
—U ee - 8 sin(y) + sin(2-))- dP C27 
dP 4Eb] 
and 
d 3-P*.a RD? . 
= cee ~ 8 sin(jz) + sin(2-yx))-dR C.28 
dR ¢ El] 





Summing equations C.25 - C.28 and simplifying, the incremental strain energy in the 


bourdon due to bending may be expressed as 

















Paw R?D 
des ep - 8 sm(u) + sin(2-)) x 
1 
4E°I 


nee r} aR e | +R-D-dP +15(P-D-dR) 


da db 








P-Rb(E-K) 
Ea’ -K-b’ 


29 


where 


qd pet C.30 


D 
da 


a’.(K - E)-(3-(a’ - 67) - 2-2?) 


= ~(%E-K) 
(a -b) 








and 


o pet 


db? 


2-a-b-(K -E) 
sare ange ak -~ 
(a’ - b’) 








as a result of taking the partial derivatives of equation C.21 with respect to the semi- 


major and minor axes, respectively. 
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APPENDIX D. DEVELOPING EQUATIONS FOR TIP TRAVEL 


Two equations are used to develop the equation for Bourdon tip displacement. 
The first equation is based on the assumption made in this analytical study that the length 
of the elastic element remains constant during pressurization. Based on this assumption, 


the following is true: 


R= Constant D.1 


and 


Be dR +R. du=0 D2 
where incremental tip travel, dR, is expressed as 


-R-d. 
ee D3 
be 


Referring to Figure 6, the variables R and us are measurable quantities but du must 
be calculated. One way to do so is to apply Castigliano’s Second Theorem as defined in 
Ref. 11 which states that 

L 


M 
a= —_—: do 
EI dM, 


dx D.4 





cal 


where M, is an imaginary moment and E’ is Young’s modulus. This equation is helpful 
because, in beam theory, @ represents the change in the slope of the elastic member at 


the point of an applied imaginary moment. Geometry then dictates that 
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sae D5 


thereby establishing the second equation necessary to calculate tip deflection. Modeling 
the Bourdon as a beam, this relationship makes sense because as the radius of curvature 
at the tip of the Bourdon becomes larger with pressure, the change in slope becomes 
larger and the angle 4 describing the length of the elastic element becomes smaller. 
Therefore, once ¢, is known, bourdon tip displacement is calculated using equation D.3. 
Continuing with the application of Castigliano’s second theorem, a dummy 
moment M/, is applied in the direction chosen positive for @ so that the bending moment 


in equation D.4 becomes 


=~ 


=-P-aR-D-(1 - cos@)-k- M-k D6 
and 
— Me-1 D.7 


Making the appropriate substitutions into equation D.4, setting 14, equal to zero, and 


integrating, 


Bh 
P-aR- -aR- 
a aRD ee a aRD 


- ToT ~ sine) 
EVI -6 Evi 








i 


D.8 
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A positive value of @ means the tip angle of rotation has the same sense of direction as 
the applied dummy moment. Because @ is a function of the dimensions a and 6 (a 
function of 5 because D is a function of 5), an iterative solution for tip travel incorporates 


the cross-sectional deformation that occurs simultaneously with straightening. 
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APPENDIX E. FORTRAN SOURCE CODE 


Program Bourdon 


This program calculates tip deflection in low pressure, C-type Bourdon 
elements for user-defined properties, geometries, and pressures. 


Definition of variables: 


a = semi-major axis of the elliptical thin-walled tube (ID) 
b = semi-minor axis of the elliptical thin-walled tube (ID) 
t = wall thickness of elliptical tube 

R = tip radius of curvature of the elliptical tube 

p = internal pressure applied to the tube 


mixs = J, the area moment of inertia of the tube for cross-sectional 


deformation 
mib = I, the area moment of inertia of the tube for bending 
mu = arclength of the elliptical tube 
epr = elastic modulus of the tube material 
nu = poisson's ratio of the tube material 
rad = input variable for the initial radius of curvature 
rin = initial radius of curvature expressed in meters 
da = incremental change in a 
db = incremental change in b 
dR = incremental change in R 
dp = incremental internal pressure applied to the tube 
dmu = incremental change in mu 


smk = k, defined as the square root of 1-(b/a)’2 
bk = K, Jacobian Elliptic Integral of the 1st kind 
be = E, Jacobian Elliptic Integral of the 2nd kind 


cl =a constant in the expression for capc 

c2dn = the product of a constant and the elliptic integral of 
dn(u)e-4 in the expression for cape 

c3dn = the product of a constant and the elliptic integral of 
dn(u)e-6 in the expression for cape 


c4dn = the product of a constant and the elliptic integral of 
dn(u)e-8 in the expression for capc 
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D =a constant in the net moment equation expressed as a 
function of a, b, K, and E 

dDda = partial derivative of D wrt a 

dDdb = partial derivative of D wrt b 


capa = a constant expressed as a function of parameters a, b, 
P, and R used in Subroutine Bndenrg 

capb = a constant expressed as a function of parameters a, b, 
P, and R used in Subroutine Bndenrg 

capc = a constant expressed as a function of dimensions a and 
b used in Subroutine Xsecenrg 

capj = an expression representative of incremental external work 

capl = the ratio of capj to capa used in the expression for capy 

capx = the ratio of capc to capa used in the final expression 
for db 

capy = an expression representing the difference in capb and 
capl used in the final expression for db 

capz = a constant expressed as a function of parameters a, b, 
P, and R used in Subroutine Bndenrg 

thetai = bourdon tip rotation due to applied dummy moment 


dV = volumetric change in the tube, dV = vol*db 
duxs = change in internal energy due to a change in the 
cross-section of the elliptical tube as a result 
of an applied internal pressure 
dub = change in internal energy due to bending of the 
elliptical tube 


* €£ & FE Fe EF HE HH KK HE HF HK HK HF HF HE KH HF KE RK HE HK HH HH HF EF KH K 


implicit real*8 (a-h,o-z) 
real*8 mixs, mib, mu, nu, maj, min 


common a,b,smk,bk,be,epr,R,mixs,mib,mu,nu,p,capc,capa,capb,capz,cl, 
+ ¢2dn,c3dn,c4dn,D,dDda,da,dDdb,pi,dp,dR,db,dmu,thetai,vol,capj, 
+ capl,capx,capy 


open (7,file='data',status='old') 
read (7,*)maj,min,wl,rad,mu,dp,pmax,epr,nu 
close (7) 


* 


* Dimensions converted to meters. 
* 


a = maj/1000.0 
b = min/1000.0 
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* %* &€ & F 


* &* & & 


t = wi/1000.0 
rin = rad/1000.0 
R=nin 


call capke 
pi = 3.141592654 


vol = pi*a*(((a**2 + b**2)*be - 2.0*bk*b**2)/(be*a**2 
+ - bk*b**2)) 


smk = sqrt(1.0 - (b/a)**2) 
mi = (pi/4.0)*((a + t)*(b + t)**3 - a*b**3) 


This loop increases the internal pressure to a pre-determined 
maximum pressure in increments of dp so that associated changes 
in the variables a, b, R, and mu may be considered in the problem. 


Open(4, file=data',status='new’) 


Do 10 p=0.0,pmax - dp,dp 
call capke 


vol = pi*a*(((a**2 + b**2)*be - 2.0*bk*b**2)/(be*a**2 
+ - bk*b**2)) 


call xsecenrg 
call dercapd 
call dummy 
call tip 

call bndenrg 


capj = (p + dp)*vol 
capl = capj/capa 
capx = capc/capa 
capy = capb - capl 


The final expression for db is written as a quadratic equation for 
both roots of the equation. 


db = (-capy + sqrt(capy**2 - 4.0*capx*capz))/(2.0*capx) 


db = (-capy - sqrt(capb**2 - 4.0*capx*capz))/(2.0*capx) 
da = (a*b*(be - bk)/(be*a**2 - bk*b**2))*db 
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* 
* 





Variables incremented through the loop. 


a=at+da 

b=b+db 

R=R+dR 

mu =mu+ dmu 

smk = sqrt(1.0 - (b/a)**2) 

mixs = t**3/12.0 

mi b= (pi/4.0)*((a + t)*(b + t)**3 - a*b**3) 


10 continue 


* € * & & 


close (4) 


stop 
end 


Subroutine Capke 


This subroutine calculates K and E, elliptic integrals of the first 
and second kind, respectively. K and E are functions of k which is 
itself a function of the major (2a) and minor (2b) axes of the 
ellipse. 


implicit real*8 (a-h,o-z) 
real*8 mixs, mib, mu, nu 


common a,b,smk,bk,be,epr,R,mixs,mib,mu,nu,p,capc,capa,capb,capz,c 1, 
+ ¢2dn,c3dn,c4dn,D,dDda,da,dDdb,pi,dp,dR,db,dmu,thetai,vol,capj, 
+ capl,capx,capy 


ang = asin(smk) 

eps = 1.0e-8 

a0 = 1.0 

bO = cos(ang) 

c0 = sin(ang) 

rmult = 1.0 

stor = cO**2 

bk = 1.5707963 

be = bk*(1.0 - 0.5*stor) 
89 if (abs(cO).It.eps) goto 99 

al = 0.5*(a0 + b0) 


80 








* ¥ &¥ & FE KF He F 


* © & 





b1 = sqrt(a0*b0) 

cl = 0.5*(a0 - b0) 

rmult = 2.0*rmult 

stor = stor + rmult*(c1**2) 
bk = 1.5707963/al 

be = bk*(1.0 - 0.5*stor) 

a0 =al 

b0 = bi 

c0=cl 

goto 89 


99 _ continue 


return 
end 


Subroutine Xsecenrg 


This subroutine calculates the change in internal energy due to a 
change in the cross-sectional area of an elliptic tube when an 
internal pressure (greater than external or ambient pressure) is 
applied to the tube. The change in internal energy occurs when 

the elliptical cross-section becomes more cicular with applied pres- 
sure. 


implicit real*8 (a-h,o-z) 
real*8 mixs, mib, mu, nu 


common a,b,smk,bk,be,epr,R,mixs,mib,mu,nu,p,capc,capa,capb,capz,c1, 
+ ¢2dn,c3dn,c4dn,D,dDda,da,dDdb,pi,dp,dR,db,dmu, thetai,vol,capj, 
+ capl,capx,capy 

call jacconst 


capc = (2.0*a*epr*mixs*c1*(c2dn + c3dn + c4dn))/(1.0 - nu**2) 


The change in internal energy due to a change in the cross-section 
is defined as duxs = capc*db**2 


return 
end 
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ee ee & 





Subroutine Jacconst 


This subroutine simplifies the expression needed to calculate the 
change in internal energy due to a change in the tube's cross- 
section. C1 is a constant for specific major and minor axes of 
the elliptical tube. Constants c2nd, c3dn, and c4dn are 

also specific for given tube axes but also include the evaluation 
of the elliptic trigonometric function dn(u) raised to the -4, -6, 
and -8 powers, respectively. 


implicit real*8 (a-h,o-z) 
real*8 mixs, mib, mu, nu 


common a,b,smk,bk,be,epr,R,mixs,mib,mu,nu,p,capc,capa,capb,capz,cl, 
+ ¢2dn,c3dn,c4dn,D,dDda,da,dDdb,pi,dp,dR,db,dmu,thetai,vol,capj, 
+ capl,capx,capy 


cl = (b/a)**4*(1.0/(be*a**2 - bk*b**2)**2) 
c2 = (((a**2 + b**2)/b**2)*be + bk)**2 

c3 = -6.0*be*(((a**2 + b**2)/b**2)*be + bk) 
c4 = 9.0*be**2 


dn2 = (2.0*be*(2.0 - smk**2)/(1.0 - smk**2) - bk)/(3.0*(1.0 - 
a smk**2)) 
dn3 = (4.0*(2.0 - smk**2)/(15.0*(1.0 - smk**2)**2))*(2.0*be*(2.0 
“+ - smk**2)/(1.0 - smk**2) - bk) - (3.0*be/(5.0*(1.0 - smk 
ie #*2)**2)) 
dn4 = (24.0*(2.0 - smk**2)**2/(105.0*(1.0 - smk**2)**3) - 5.0/ 
na (21.0*(1.0 - smk**2)**2))*(2.0*be*(2.0 - smk**2)/(1.0 - 
+ smk**2) - bk) - (18.0*be*(2.0 - smk**2)/(35.0*(1.0 - smk**2)**3)) 
c2dn = c2*dn2 
c3dn = c3*dn3 
c4dn = c4*dn4 


return 
end 


Subroutine Dummy 


This subroutine calculates the change in the slope at the free end 

of the tube using Castigliano's second theorem. The change in the 
slope is determined by applying a dummy moment at the tube tip and 
incorporating that moment into the internal energy equation for a 
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* beam. Thetai is found by integrating the internal energy equation and 
then setting the dummy moment equal to zero. 


* & &€ 


* * &€ & & * 


implicit real*8 (a-h,o-z) 
real *8 mixs, mib,nu, mu 


common a,b,smk,bk,be,epr,R,mixs,mib,mu,nu,p,capc,capa,capb,capz,c1, 
+ ¢2dn,c3dn,c4dn,D,dDda,da,dDdb,pi,dp,dR,db,dmu,thetai,vol,capj, 
+ capl,capx,capy 


thetai = ((p + dp)*a*(R**2)*D*(mu - sin(mu)))/(epr*mib) 


Thetai should be a positive value to indicate that the angle of 
rotation has the same sense of direction as the applied dummy moment. 


return 

end 

Subroutine Tip 
Once thetai is known, the tip deflection of the tube can be 
calculated by assuming that mu*R = constant. In doing so, 
mu*dR + R*dmu = 0 must be true. Geometry shows that mu has the 


same magnitude but opposite sense of thetai. 


implicit real*8 (a-h,o-z) 
real *8 mixs, mib, nu, mu 


common a,b,smk,bk,be,epr,R,mixs,mib,mu,nu,p,capc,capa,capb,capz,c1, 
+ ¢2dn,c3dn,c4dn,D,dDda,da,dDdb,pi,dp,dR,db,dmu,thetai,vol,capj, 
+ capl,capx,capy 

dmu = -thetai 

dR = (-R*dmu)/mu 


return 
end 
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Subroutine Bndenrg 


This subroutine simplies the expression needed to calculate the 
change in internal energy due to bending in the tube. Capa, capb, 
and capz are all constant for the given parameters but change as the 
parameters change. 


implicit real*8 (a-h,o-z) 
real *8 mixs, mib, nu, mu 


common a,b,smk,bk,be,epr,R,mixs, mib,mu,nu,p,capc,capa,capb,capz,cl, 
+ c2dn,c3dn,c4dn,D,dDda,da,dDdb,pi,dp,dR,db,dmu, thetai, vol,capj, 
+ capl,capx,capy 


capa = (((p + p)*(a**2)*(R**2)*D)*(6.0*musin(mu)+sin(2.0*mu)))/ 
i (4.0*epr*mi) 

capb = (((p + dp)*R*b*(be - bk)/(be*a**2 - bk*b**2))*(a*dDda + 
+ — D))+ ((p + dp)*R*dDdb) 

capz = (R*D*dp) + (1.5*(p + dp)*D*dR) 


The change in internal energy due to bending is defined as 
dub = capa*(capb*db + capz) 


return 
end 


Subroutine Dercapd 


This subroutine simplifies the expression for the net moment acting 
on the tube as a result of the applied internal pressure. D is the 
simplifying term in the moment equation. dDda and dDdb are the 
partial derivatives of D with respect to a and b, respectively. 


implicit real*8 (a-h,o-z) 
real*8 mixs, mib, mu, nu 


common a,b,smk,bk,be,epr,R,mixs,mib,mu,nu,p,capc,capa,capb,capz,cl, 


+ ¢2dn,c3dn,c4dn,D,dDda,da,dDdb,pi,dp,dR,db,dmu, thetai,vol,capj, 
+ capl,capx,capy 
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D = a*((4.0/3.0)*((bk - be)/(smk**2) + (2.0*be - bk) - pi*(b/ 
a)) 

dDda = (4.0/3.0)*(((a**2*(bk - be)*(3.0*(a**2 - b**2) - 2.0*a** 
2))(a**2 - b**2)**2) - (2.0*be - bk)) 


dDdb = (4.0/3.0)*((2.0*a**3*b*(bk - be)/(a**2 - b**2)**2) - pi) 


return 
end 


85 











86 

















LIST OF REFERENCES 


1. Hartland, P.W., Pressure Gauge Handbook, Marcel Dekker, Inc., 1985. 
2. Heise, Dresser Industries Instrument Division, Newtown, CT. 
3. ASME B.40.1, Gauges - Pressure Indicating Dial Type - Elastic Element, 1991. 


4. Jimenez-dominguez, H., Figueroa-Lara, F., and Galindo, S., “Bourdon Gauge 
Absolute Manometer”, /nstituto Nacional de Investigaciones Nucleares, pp. 499, 


Nov., 1985. 
5. Stewart Warner Instrument Corporation, Mt. Prospect, IL. 
6. Wika Instrument Corporation, Lawrenceville, GA. 


7. Giacobbe, J. B. and Bounds, A. M., “Material Selection Factor Significant in 
Bourdon Tubes”, Journal of Metals, pp. 1147, Nov., 1952. 


8. Wolf, Alfred, “An Elementary Theory of the Bourdon Gage”, Journal of Applied 
Mechanics, pp. A-207, Sep., 1946. 


9. Bowman, F., Introduction to Elliptic Functions with Applications, English 
Universities Press LTD, London, 1953. 


10. Byrd, P.F. and Friedman, M.D., Handbook of Elliptic Integrals for Engineers and 
Scientists, 2nd Ed., Springer-Verlag Berlin, 1971. 


11. Logan, D.L., Mechanics of Materials, Harper Collins Publishers, 1991. 


12. Andreeva, L.E., Elastic Elements of Instruments, Daniel Davey & Co., Inc., 
N.Y., 1966. 














13. Ribreau, C., Naili, S., Bonis, M., and Langlet, A., “Collapse of Thin-Walled 
Elliptical Tubes for High Values of Major-to-Minor Axis Ratio”, Journal of 
Biomechanical Engineering, pp. 432, Vol. 115, Nov., 1993. 


14. Gradshtayn, I'S. and Ryzhik, I.M., Table of Integrals, Series, and Products, 
Academic Press, Inc., San Diego, CA, 1980. 


88 








INITIAL DISTRIBUTION LIST 


No. Copies 
1. Defense Technical Information Center Z 
8725 John J. Kingman Rd., STE 0944 
Ft. Belvoir, VA 22060-6218 


2. Library, Code 13 2 
Naval Postgraduate School 
Monterey, CA 93943-5101 


3. Department Chairman, Code ME 2 
Department of Mechanical Engineering 
Naval Postgraduate School 
Monterey, CA 93942-5000 


4. Professor Ranjan Mukherjee 1 
Department of Mechanical Engineering 
Naval Postgraduate School 
Monterey, CA 93942-5000 


5. Curricular Officer, Code 34 1 
Department of Mechanical Engineering 
Naval Postgraduate School 
Monterey, CA 93942-5000 


6. LT Cynthia D. Conway 1 


33 McCauley Rd. 
Travelers Rest, SC 29690 


89 





